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journal  articles.  The  authors  wish  to  thank  Lt  Col  Ben  A.  Loving, 

Capt  S.  A.  Shackelford,  Maj  Joel  W.  Beckmann,  Capt  Charles  L.  Hussey, 
and  Lt  Larry  Wright  of  AFWL  for  helpful  discussions.  In  addition,  thanks 
are  due  to  D.  G.  Reynolds  and  B.  J.  Darcy  for  typing  the  manuscript. 
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i.  introduction 


Many  chemical  problems  are  amenable  to  theoretical  treatment  with 
rrolecular  orbital  techniques.  These  calculations  have  new  evolved  to 
a point  where  they  can  be  reliably  used  to  predict  molecular  geometries, 
ionization  potentials,  potential  energy  curves,  heats  of  formation,  and 
reaction  activation  energies.  The  goal  of  this  work  unit  is  to  apply 
available  molecular  orbital  programs  to  chemical  problems  of  Air  Force 

interest. 

Several  of  these  problems  have  been  investigated.  First,  applica- 
tion of  available  semiempirical  programs  to  calculation  of  potential 
curves  of  diatomic  molecules  of  laser  interest  was  attempted.  The  most 
advanced  of  the  smiempirical  programs,  MNDO,  developed  by  Dewar  and 
coworkers  at  the  University  of  Texas1,  was  employed  to  calculate  poten- 
tial curves  for  the  molecule  NF.  This  report  concludes  the  effort  to 
use  semiempirical  programs  to  do  calculations  on  diatcmics  of  laser 

interest. 

The  Air  Force  Rocket  Propulsion  Laboratory  needs  accurate  molecular 
geometries  and  heats  of  formation  for  boron-  and  aluminum-containing 
species2.  These  heats  of  formation  are  used  as  inputs  to  rocket  engine 
performance  codes.  MINDO/3  has  been  used  to  calculate  heats  of  formation 
and  geometries  of  boron-containing  species.  Vfork  is  currently  underway 
to  use  MNDO  to  calculate  heats  of  formation  and  geometries  of  boron-  and 
aluminum-containing  species. 

Experimental  studies  on  the  decomposition  rates  and  mechanisms  of 
explosives  have  been  underway  in  the  laboratory  for  quite  seme  time  ' . 


MINDO/3  and  now  MNDO  are  being  employed  to  calculate  energies  of  sus- 
pected decomposition  products  from  TOT.  These  energies  can  be  used  to 
estimate  activation  energies  for  possible'  decomposition  mechanisms. 

Thus  the  calculations  can  be  used  as  an  indication  of  mechanism.  The 
eventual  goal  for  this  work  is  the  correlation  of  explosives  properties 
with  the  calculated  molecular  orbital  properties. 

It  would  be  of  great  utility  to  be  able  to  predict  oxidation  and 
reduction  potentials  for  species  in  electrochemical  systems.  Since 
oxidation  potentials  are  related  to  ionization  potentials  which  can  be 
calculated  via  Koopmans'  theoran^  from  molecular  orbital  energies,  one 
would  expect  to  be  able  to  correlate  MO  energies  with  oxidation  potentials. 
The  same  relationship  holds  between  reduction  potentials  and  electron 
affinities.  Attempts  to  make  these  correlations  have  been  reasonably 
successful  with  Huckel  theory^.  In  this  effort,  MNDO  was  anployed  to 
predict  oxidation  and  reduction  potentials  for  several  aromatic  molecules. 
In  addition,  electron  affinities  for  several  aluminum-containing  species 
predicted  to  be  present  in  a NaCl-AlCl^  molten  salt  solution  have  been 
calculated  to  determine  which  species  are  more  acidic. 


2 


Several  semiempirical  methods  have  been  used  in  this  work.  All  are 
based  on  a Hartree  ?ock  self-consistent  field  technique  including  only 
the  valence  electrons.  Both  open-shell  and  closed-shell  molecules  have 


been  calculated. 

The  rrost  basic  program  used  was  Quantum  Chemistry  Program  Exchange 
(QCPE)  Program  Number  261,  Complete  Neglect  of  Differential  Overlap 
Method  Including  Third  to  Elements.  The  program  is  an  extension  of  one 
developed  by  Pople  and  Beveridge7.  It  calculates  molecular  energies, 
dipole  noments,  nolecular  orbitals,  and  molecular  orbital  energies.  It 
uses  a fixed  molecular  geometry  as  input  and  is  not  capable  of  optimizing 
geometries.  It  has  been  used  primarily  in  this  work  in  its  MX)  ( Inter- 
mediate Neglect  of  Differential  Overlap)  option  to  calculate  ESR 

splitting  constants  for  radical  systems. 

Ihe  CNDO/S  program  developed  by  Jaffe  and  Del  Bene  was  developed 
for  calculating  electronic  transition  energies.  It  is  based  on  the  CNDO 
method  but  contains  parameters  specifically  to  give  a good  fit  for  elec- 
tronic transition  energies  of  aromatic  molecules.  It  has  been  used  in 
this  work  to  calculate  transition  energies  for  TOT.  It  is  only  capable 

of  calculating  closed  shell  molecules. 

Ihe  MINDO  series  of  programs  developed  by  Dewar  and  coworkers  at  the 
University  of  Itexas9 ' 10 ' 11  have  developed  into  useful  analytical  tools 
for  studying  reaction  paths,  reaction  energies,  and  molecular  geometries. 
Ihe  MINDO/3  version11  has  been  used  in  the  early  stages  of  this  work. 
MINDO/3  has  been  shown  to  give  reliable  geometries  and  heats  of  formation 
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for  a wide  variety  of  molecules  ‘containing  carbon,  oxygen,  nitrogen, 
and  hydrogen  ' . Parameters  for  fluorine,  chlorine,  silicon,  phos- 

phorus, and  sulfur  have  also  been  included  in  order  to  do  calculation 

14  15 

on  molecules  containing  these  atans  ' . Recently  a more  sophisticated 

semiempirical  program  called  MNDO  based  on  an  NDDO  (neglect  of  diatonic 
differential  overlap)  scheme  has  been  developed  by  Dewar's  group'*'.  It 
has  been  used  exclusively  for  the  molecular  energy  calculations  since  it 
was  received  in  the  spring  of  1977.  It  is  new  capable  of  doing  calcu- 
lations on  molecules  containing  C,N,0,H,F,Cl,Si,P,S  and  Al.  Results 
with  MNDO  have  been  shown  to  be  better  than  those  from  MINDO/3. 


i . 

/>  ' 

’ - 
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III.  RESULTS  AND  DISCUSSION 


Potential  Curves  of  Diatomic  Molecules 

Accurate  potential  energy  curves  for  electronic  states  of  diatonic 
molecules  are  necessary  for  a consideration  of  these  molecules  as  lasing 
species.  In  addition,  the  potential  curves  are  used  as  input  to  calcu- 
lations of  the  reaction  rates  between  molecules.  Ab  Initio  programs 
have  been  used  to  calculate  these  curves  when  experimental  data  has  not 
been  sufficient  to  construct  RKR  potential  curves.  If  available  semi- 
empirical  programs  could  be  used  instead  of  the  ab  initio  programs, 
calculation  of  these  potential  curves  would  be  much  easier  and  quicker. 

As  a preliminary  screening  tool,  we  used  the  semiempirical  programs 
previously  described  to  calculate  the  ground  state  curve  for  HF.  The 
CNDO/S  program  gave  no  minimum  in  the  curve;  this  failure  is  not  surpris- 
ing because  of  its  parameterization  for  transition  energies.  The  CNDO/2 
program  (both  in  its  CNDO  and  INDO  options)  and  the  MINDO/3  program  gave 
curves  of  the  appropriate  Morse  shape,  but  the  results  were  not  good  at 

O 

intemuclear  distance  greater  than  2A.  All  of  these  correctly  predicted 

O 

the  equilibrium  intemuclear  distance  to  be  slightly  less  than  one  A. 
Since  the  CNDO/2  program  and  MINDO/3  are  not  capable  of  doing  calcula- 
tions on  excited  states,  MNDO  was  chosen  to  do  calculations  on  all 
possible  excited  state  configurations  of  the  NF  molecule. 

The  NF  molecule  has  8 valence  orbitals  (four  from  each  atom)  and  12 
valence  electrons.  The  valence  shell  configuration  is 


2 2 2 2 2 1 1 
°1  °2  °3  irl,-lTrl,+lir  2,-l1T  2,+l 
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For  the  tt  orbitals,  the  first  subscript  referes  to  the  unperturbed  energy 
and  the  second  refers  to  the  value  of  the  magnetic  quantum  number.  Thus, 


tt^  and  are  a degenerate  set  of  tt orbitals;  likewise  for 

and  ^2+1  <I^e  MND0  program  allowed  16  configuration  interaction  states. 

These  are  listed  below; 

(1)  Oj  o2  03  n1/_1  *li+1  ”2,-l  ^2 ,+l  °4° 

(2)  a12<,22o32,li.12,li+12,2>.12,2i+1o4- 

2 2 2 2 

(3)  a1  a2  a3  it1>+1  *2,-1  'n2,+l  a4 

2 2 2 2 1 1 2 

(4)  al  °2  °3  *1,-1  *1,+1  1T2,-1  *2,+l  °4 

(5)  cr12a22a  3 27T1>_117:1>+12TT2f_11Tr2f+12c4° 

(6)  ox  c2  a3  _x  Trlf+1  "2,-1  ^2,+l  °4° 

(1)  a12a21a327Tlf_127rlf+12Tr2f_11TT2/+12a4° 

(8)  CT1  c2  a3  TT1^_1  7r1/+1  ^2,-1  *2,+l  °4° 

(9)  a2  a3  tt-l  _x  *1/+1  ^2,-l  1T2,+l°a4 

(10)  a1  a2  a3  tTj^  *lf+1  ff2,-l  *2,+l  a4° 

(11)  o2  a3  •'T1^_1  ^1;+1  ir2/-l  1T2,+1  °4° 

(12)  ox  a2  o3  ^1/+1  1t2,-1  *2,+l  °4° 

(13)  a1  a2  a3  *lf+1  *2,-1  *2,+l  a4° 

(14)  ax  a2  a3  ir2,-l  *2,+l  a 4° 

(15)  a2  o3  if1#+1  1T2,-1  *2,+l  0 4° 

222  2 2 0 °2 

(16)  Oj.  a2  a3  ff1>+1  *2,-1  n2,+l  0 4 


All  sixteen  of  these  configurations  were  done  in  the  singlet  state 
manifold.  For  the  triplet  state  manifold,  configurations  must  have 
unpaired  electrons.  Configurations  (1),  (4),  (11),  and  (16)  can  be  mixed 
to  form  I+  states.  Configurations  (2),  (3),  (5),  and  (10)  are  of  A 
synmetry,  and  the  remainder  are  of  it  synmetry.  Thus,  from  these  16 
configurations  we  expect  4E+,  2A,  and  4tt  states  in  the  singlet  manifold. 

In  the  triplet  manifold , we  expect  3 E+,  1A,  and  4tt  states  frcm  the  13 
configurations.  The  problem  with  identification  of  these  states  frcm 
the  MNDO  output  is  that  the  Hamiltonian  is  not  synrne try-adapted , and 
therefore  states  of  different  synmetry  do  mix.  For  some  situations,  the 
mixing  is  not  substantial  enough  to  preclude  state  identification;  in 
other  cases,  it  is. 

The  program  did  predict  the  triplet  state  to  be  the  most  stable,  and 

O 

the  equilibrium  intemuclear  distance  was  predicted  to  be  1.220  A.  The 
heat  of  formation  was  predicted  to  be  31.6  kcal/mole.  No  potential  curves 
will  be  presented  because  of  the  problems  in  identifying  and  separating 
states  of  different  synmetry. 


Rocket  Plume  Species 


MINDO/3  was  used  to  make  calculations  on  a number  of  boron-aontain- 
ing  species  in  order  to  determine  its  accuracy  for  these  molecules. 

Both  ground-state  geometries  and  heats  of  formation  were  predicted.  The 
results  are  given  in  Table  I. 

Table  I shows  that  the  worst  predictions  are  made  for  the  diatonic 
species.  Both  the  heat  of  formation  and  the  intemuclear  distance  are 
considerably  off  for  both  BF  and  BO.  Calculation  of  the  heat  of  forma- 
tion using  experimental  intemuclear  distances  only  improved  the  results 
slightly.  The  predictions,  particularly  the  heats  of  formation,  get 
considerably  better  for  the  larger  molecules.  In  the  case  of  BF^,  the 
MINDO/3  results  are  excellent,  because  BF^  was  probably  one  of  the 
molecules  used  to  deduce  the  semiempirical  parameters  for  boron. 

The  geometry  of  BF2  is  not  well  known.  The  bond  distance  is  an 
assumed  one,  and  the  112°  bond  angle  is  deduced  from  ESR  studies. 

MINDO/3  predicts  a linear  molecule,  and  its  prediction  for  the  heat  of 
formation  is  considerably  off.  When  the  experimentally-deduced  geometry 
is  used,  the  MIND0/3-predicted  heat  of  formation  improves  tremendously; 
it  becomes  -119.2  kcal/mole,  in  much  better  agreement  than  the  previous 
result.  This  better  agreement  is  a good  indication  that  BF2  is  indeed 
a bent  and  not  a linear  molecule. 

The  predictions  for  B02  and  FBO  are  both  excellent.  The  assumed 
geometry  for  FBO  agrees  well  with  that  predicted  by  MINDO/3.  Both 
heats  of  formation  are  almost  within  the  experimental  error. 

EF20  caused  some  problems  for  the  MINDO/3  calculation.  When  the 
geometry  was  optimized,  the  optimization  resulted  in  a molecule  with  two 
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TABLE  I 


MINDO/3  vs.  Experimental  Parameters 


for  Boron-Cont 


Molecules 


MOLECULE 

PARAMETER 

EXPERIMENTAL1 

MINDO/3 

BF 

B-F  Bond  Distance 
AHf,  25°C 

1.265A 

-27.7+3.3  kcal/mole 

1.196A 

-107.5  kcal/mole 

BO 

B-0  Bond  Distance 
AHf,  25°C 

1.2049A 

0+2  kcal/mole 

1.145A 

-55.9  kcal/mole 

^2 

B-F  Bond  Distance 
F-B-F  Bond  Angle 
AHf,  25°C 

° 2 

(1.30A) 

112° 

-141.0+3.0  kcal/mole 

1.264A 

180.0° 

-207.3  kcal/mole 

bo2 

B-0  Bond  Distance 
O-B-O  Bond  Angle 
AHf,  25°C 

1.263A 

180° 

-68+2  kcal/mole 

1.206A 

179.9° 

-72.0  kcal/mole 

FBO 

F-B  Bond  Distance 
B-O  Bond  Distance 
F-B-O  Bond  Angle 
AHf/  25°C 

(1.30|) 

(1.20A) 

(180°) 

-144+3  kcal/mole 

1.304A  0 

1.1.184A 

180.0° 

-145.4  kcal/mole 

bf3 

B-F  Bond  Distance 
F-B-F  Bond  Angle 
AHff  25°C 

1.307+.002A 

120°  ~ 

-271.42+.4  kcal/mole 

1.302A 

120.0° 

-270.2  kcal/mole 

BF2° 

BF  Bond  Distance 

BO  Bond  Distance 
O-B-F  Bond  Angle 
F-B-F  Bond  Angle 
AH-,  25°C 

1.30+.05i 

1.40+.05A 

117+5° 

126+5° 

-200+15  kcal/mole 

1.318| 

1.247A 

121.5° 

117.0° 

-189.4  kcal/mole 

il 


"'"All  data  taken  from  JANNAF  Thermochemical  Data  Tables. 

^Brackets  around  a parameter  mean  that  it  is  an  assumed  parameter. 
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different  O-B-F  bond  angles.  The  results  given  are  for  those  of  a 
synmetry— enforced  calculation  which  forced  both  O-B-F  bond  angles  to 
the  same  value.  In  either  case,  the  predicted  geometry  was  considerably 
r arrived  from  the  experimental  geometry.  However,  in  both  cases  the 
predicted  heats  of  formation  was  within  the  experimental  error  of  the 
experimental  value. 

Calculations  are  now  continuing  with  MNDO  on  this  same  set  of  mole- 
cules. In  addition,  aluminum-containing  species  will  also  be  studied. 
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ELECTROCHEMICAL  STUDIES 


r 

■I-, 


A series  of  aromatic  molecules  have  been  calculated  with  MNDO  in 
order  to  correlate  the  observed  electrochemical  oxidation  potentials 
with  calculated  highest  occupied  molecular  orbital  (HOMO)  energies. 
These  oxidation  potentials  were  measured  in  acetonitrile/0.5  sodium 
perchlorate  solution  vs.  the  Ag/O.IM  Ag+  electrode16.  The  correlation 
is  shown  in  Figure  1. 

The  most  obvious  feature  of  Figure  1 is  the  excellent  correlation 
between  the  oxidation  potential  and  the  HOMO  energy  for  this  series  of 
compounds.  A least  squares  fit  (shewn  on  the  figure)  gives  a slope  of 
-1.08  and  an  intercept  of  -7.17.  Biphenyl  was  not  included  in  the 
least  squares  fit  because  of  a geometry  optimization  problem.  When 
biphenyl  is  completely  optimized,  the  angle  between  the  rings  is  pre- 
dicted to  be  12°.  The  experimental  gas  phase  angle  is  about  50°.  Both 
results  are  shown  on  the  figure.  The  results  with  the  50°  angle  are 
better  to  such  a degree  that  one  is  tempted  to  use  this  correlation  as 
a predictive  tool  to  deduce  the  correct  ring  angle. 

These  results  are  very  encouraging,  not  only  from  the  standpoint  of 
the  straight  line,  but  also  because  of  the  slope  being  very  close  to  a 
value  of  1.00.  This  is  a strong  indication  that  for  these  systems  the 
solvent  interactions  are  doing  very  little  to  perturb  the  molecular 
orbital  energies,  or  at  least  they  are  all  being  perturbed  in  a like 
manner.  This  correlation  will  be  used  to  predict  oxidation  potentials 
for  anthracene  derivatives.  Deviations  frem  the  correlation  will  be 
used  to  indicate  possible  extensive  solvent-solute  interactions. 
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MNDO  calculations  are  continuing  on  aluminum  chloride  species  of  an 
AlCl3/NaCl  molten  salt  system.  The  principal  aim  is  to  calculate 
electron  affinities  of  these  species  to  determine  which  are  more  likely 


to  remove  electrons  from  aromatic  molecules. 


HOMO  ENERGY  (eV) 


-8.0 


-8.5 


BIPHENYL  (50°) — ► 4 


-9.0 


BIPHENYL  (72°)  — ► * 


-9.5  4 ■■  | V 

•75  1.0  1.25  1.5  1.75  2.0 

OXIDATION  POTENTIAL  (V) 


FIGURE  1.  CORRELATION  BETWEEN  OXIDATION 
POTENTIALS  AND  HOMO  ENERGIES. 
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TNT  THERMAL  DECOMPOSITION  MECHANISM 


Both  MINDO/3  and  MNDO  have  been  used  to  calculate  TOT  and  suspected 
decomposition  products  in  order  to  deduce  the  mechanism  for  TOT  decompo- 
sition. Estimates  of  activation  energies  for  suspected  reaction  path- 
ways can  be  measured  and  compared  with  measured  experimental  values. 

TOT  decomposition  has  been  studied  experimentally  for  many  years. 

The  liquid  phase  decomposition  is  extremely  complex  and  probably  involves 
autocatalysis.  It  is  now  evident  that  the  rate-determining  step  in  the 

early  phases  of  decomposition  involves  the  breakage  of  a carbon-hydrogen 

18 

bond  of  the  methyl  group  . This  step  apparently  has  an  activation 

19 

energy  of  about  46.5  kcal/mole,  as  deduced  from  induction  period  data  . 

Thus,  the  initial  theoretical  work  has  involved  testing  of  proposed 

reaction  pathways  for  this  rate-determining  step  and  comparing  calculated 

activation  energies  for  each  pathway  with  the  measured  value. 

Quantum  mechanical  studies  of  explosives  have  been  carried  out  for 

several  years  by  Michael  Schroeder  at  the  U.S.  Army's  Ballistic  Research 
20-26 

Laboratories  . The  general  goal  of  that  work  has  been  to  correlate 
explosives  properties  with  molecular  orbital  calculational  results. 

The  program  has  been  reasonably  successful,  and  shews  the  utility  of 
these  techniques  in  correlating  explosives  properties.  The  work 
described  in  our  report  is,  to  our  knowledge,  the  first  attempt  to 
deduce  TOT  decomposition  mechanisms  from  molecular  orbital  calculations. 

It  is  generally  agreed  that  the  first  step  of  TOT  decomposition 
occurs  as  an  oxidative  attack  by  a nitro  group  on  the  methyl  group.  The 
possibility  exists  that  this  could  occur  in  either  an  inter-  or  an 
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intramolecular  fashion.  In  addition,  either  of  these  processes  could 
occur  as  either  an  insertion  step  in  which  the  oxygen  atom  inserts 
between  hydrogen  and  the  methyl  grcup  carbon  or  as  a hydrogen  atom 
transfer  from  the  methyl  group  to  the  nitro  group.  For  intramolecular 
steps,  the  former  would  result  in  a nitroso  alcohol  and  the  latter  in 
a diradical  species.  In  the  case  of  an  intermolecular  step,  the  inser- 
tion would  result  in  a nitroso  compound  and  an  alcohol,  and  the  transfer 
would  produce  two  radical  species.  The  para  NC^  group  could  also  be 
involved  in  any  of  these  steps.  These  possible  first  steps  are  shown 
in  Figure  2. 

As  a first  step  in  attempting  to  calculate  molecular  energies  along 
the  reaction  paths  mentioned  above,  nitrobenzene  was  studied^,  of 
primary  interest  was  the  correctness  of  the  geometry  prediction  in 
terms  of  the  angle  between  the  phenyl  group  and  the  nitro  group.  The 
MINDO/3  program  incorrectly  predicted  that  the  two  are  perpendicular; 
the  experimental  evidence  indicates  a planar  molecule.  MNDO  does  not 
predict  a perpendicular  conformation,  but  it  still  predicts  a non-planar 
molecule.  However,  the  energy  differences  between  these  conformations 
is  only  a few  kcal/mole.  Therefore,  any  errors  made  in  the  heats  of 
formation  for  these  nitro  group-containing  molecules  will  only  be  a few 
kcal/mole,  and  there  probably  will  be  at  least  partial  cancellation 
between  reactants  and  oroducts  of  these  errors. 

MINDO/3  was  used  to  calculate  heats  of  formation  of  TNT  and  sane  of 
the  possible  intermediates  shewn  in  Figure  2.  Species  I,  II,  III,  IV,  and 
VI  have  been  completed.  In  addition,  the  OH  radical,  hydroquinone  (VIII) , 
the  hydroquinone  radical  (IX) , and  quinone  (X)  have  been  done  in  order 
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INTERMOLECULAR  OXIDATIVE  INSERTION 


X X 


VI 


INTRAMOLECULAR  OXIDATIVE  INSERTION 


X s N02  X VII 

INTRAMOLECULAR  OXIDATIVE  TRANSFER 

FIGURE  2.  POSSIBLE  FIRST  OXIDATIVE  STEPS 
FOR  TNT  THERMAL  DECOMPOSITION. 
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to  assess  potential  reactions  between  hydroguinone  and  TOT.  The 
diradical  species  VII  was  also  attempted,  but  the  SCF  calculation  would 
not  converge.  Initially  this  was  thought  to  mean  that  it  would  not  be 
very  stable,  but  it  could  very  well  mean  that  the  initial  approximation 
to  the  density  matrix  was  not  good  enough  to  allow  convergence.  MNDO 
has  been  used  to  calculate  these  same  species  since  it  has  been  operable 
on  our  system.  To  date  the  hydroquinone  species,  OH,  TOT,  and  II  have 
been  done.  The  results  for  all  of  these  calculations  are  given  in 
Table  II.  Efforts  are  continuing  to  calculate  the  remaining  species 
with  MNDO. 

The  results  with  the  MINDO/3  method  given  in  Table  II  were  not  done 
with  complete  geometric  optimization.  All  aromatic  rings  were  assumed 
planar,  and  the  NO2  groups  ortho  to  a methyl  or  a methylene  group  were 
held  fixed  at  a 45°  angle  to  the  ring.  Para  N02  groups  were  held  fixed 
in  a coplanar  configuration  with  the  ring.  Hie  MINDO/3  method  uses 
the  half-electron  method  to  calculate  radicals. 

In  the  case  of  MNDO,  complete  geometry  optimization  was  allowed.  In 
the  case  of  the  radicals,  geometry  optimization  was  done  with  the  unre- 
stricted Hartree-Fock  procedure  (UHF) . The  final  value  for  the  heat  of 
formation  was  then  calculated  with  the  half-electron  method  utilizing 
the  UHF-optimized  geometry.  These  procedures  were  followed  because  the 
UHF  optimization  was  much  more  rapid  than  the  half-electron  optimization, 
but  the  MNDO  parameters  are  set  such  that  the  heats  of  formation  calcu- 
lated with  the  half-electron  method  are  more  accurate  than  those  calcu- 
lated with  the  UHF  method.  Sample  calculations  on  the  hydroquinone 
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Heat  of  Formation  at  25°  (kcal/irole) 


MOLECULE 

MINDO/31 

MNDO2 

OH 

16.9 

0.5 

TNT 

-11.3 

76.4 

II 

23.4 

111.5 

III 

-1.1 

- 

IV 

11.6 

- 

VII 

-84.1 

-74.0 

IX 

-49.9 

-40.7 

X 

-40.8 

-35.7 

^All  aromatic  rings  assumed  to  be  planar,  NO^  groups  ortho  to  methyl 
or  a methylene  group  held  fixed  at  a 45°  angle  to  the  ring;  para  NC^ 
fixed  coplanar  with  the  ring. 

2 

Complete  geometric  optimization  allowed,  heat  of  formation  calculated 
with  the  half-electron  method  for  radicals  using  the  UHF -optimi zed 
geometry. 
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radical  (IX)  produced  essentially  the  same  heat  of  formation  with  the 
half -electron  method  utilizing  its  own  optimized  geometry.  For  both 
MNDO  and  MINDO/3  calculations  on  quinone,  configuration  interaction  was 
included  because  of  its  possible  diradical  character. 

Table  II  gives  some  interesting  results.  For  the  quinone  species 
(VIII,  IX,  and  X),  MINDO/3  and  MNDO  results  are  reasonably  close.  How- 
ever, the  TOT  species  show  quite  a large  discrepancy  between  the  two 
methods.  It  is  interesting  to  point  out  that  the  MINDO/3  results  for 

TOT  itself  is  much  closer  to  the  true  experimental  value  of  -17.8  kcal 

27 

mole  than  the  more  advanced  MNDO  method.  The  difference  may  lie  in 
the  difference  in  the  method  of  optimization  in  each  case,  but  since 
MNDO  was  completely  optimized  geanetircally  while  MINDO/3  was  not,  one 
would  expect  MNDO  to  give  a lower  heat  of  formation  than  MINDO/3,  but 
it  does  not.  The  difference  in  energy,  however,  between  the  TNT  molecule 
and  its  derived  radicals  appears  to  be  about  the  same  for  both  methods. 
This,  the  estimated  E 's  will  be  about  the  same.  The  two  calculated 
values  for  the  OH  radical  (16.9  kcal/mole  for  MINDO/3  and  0.5  kcal/mole 
for  MNDO)  bracket  the  experimental  value  of  9.49  + .04  kcal/mole  as  taken 
from  the  JANNAF  tables. 

In  order  to  compare  with  the  measured  activation  energies,  the 
transition  state  geometry  would  have  to  be  found.  For  endothermic 
reactions,  a calculation  of  AH  for  the  reaction  will  give  a lower  bound 
to  the  activation  energy.  This  AH  for  several  potential  system  reactions 
can  be  calculated  from  Table  II.  Efforts  are  currently  underway  to 


obtain  the  transition  state  geometry  and  thus  a better  estimate  of  E^. 


One  possible  intermolecular  first  step  for  TOT  decomposition  from 
Figure  2 is: 

2TOT  II  + III 

Using  the  tabulated  heats  of  formation  of  MINDO/3  from  Table  II,  the  AH 

for  this  reaction  is  44.6  kcal/mole.  The  minimization  for  Species  III 

has  not  yet  been  completed  with  MNDO,  but  the  energy  difference  between 

TOT  and  Species  II  indicates  that  the  MNDO-calculated  for  this  step 

will  be  close  to  that  of  MINDO/3 . The  MINDO/3  value  is  certainly  close 

to  the  experimentally-measured  value  for  the  rate-determining  step  of 

4 

TOT  decomposition  of  46.5  + 1.5  kcal/mole  . 

Unfortunately,  not  all  of  the  other  species  have  been  calculated,  so 
that  estimates  of  E for  the  other  possible  first  steps  given  in  Figure  2 

d 

could  be  made.  These  calculations  are  continuing.  We  can,  however, 
estimate  an  E for  the  possible  further  degradation  of  species  III: 

cl 

III  -*  IV  + OH 

AH  for  this  process  is  calculated  to  be  29.9  kcal/mole  with  MINDO/3. 

In  addition,  AH's  for  the  intermolecular  reactions  of  the  quinone 
species  with  TOT  can  be  calculated.  The  donation  of  a hydrogen  to  TOT 
can  be  written  as: 

VIII  + TOT  -*  IX  + III 

AH  for  this  reaction  is  calculated  to  be  44.1  kcal/mole  with  MINDO/3. 

For  the  further  donation  fran  IX  to  another  TOT: 

IX  + TOT  ■+  X + III 

AH  is  calculated  to  be  19.0  kcal/mole.  The  first  donation  has  an  acti- 
vation energy  only  slightly  less  than  the  TOT-TOT  transfer,  but  only 
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r 1 

i 

needs  to  be  slightly  less  in  order  to  speed  the  decomposition  signifi- 

5 

cantly.  Also,  the  second  donation  from  the  hydroquinone  radical  species 

IX  has  a much  lower  activation  energy.  Experimentally,  the  induction 

. 

period  of  the  decomposition  decreases  markedly  when  hydroquinone  is 
added  to  TOT. 

Several  points  need  to  be  made  about  these  calucalations  used  for 
estimating  Ea  for  these  reactions.  First,  the  calculated  AH's  obtained 
thus  far  are  only  lcwer  bounds  to  the  actual  calculated  activation 
energies.  Second,  these  MO  calculations  are  on  isolated  molecules,  and 
thus  are  strictly  applicable  only  for  gas  phase  reactions.  However,  if 
intermolecular  effects  on  the  reactions  are  negligible  in  the  liquid 
phase,  then  the  calculated  energies  will  be  good  approximations  to  the 
experimental  activation  energies.  Third,  use  of  activation  energies 
alone  to  predict  relative  reaction  rates  for  different  reactions  is  not 
valid  unless  the  frequency  factors  or  preexponentials  in  the  Arhennius 
expression  are  similar. 


J 


21 


The  attempt  to  use  semiempirical  molecular  orbital  techniques  to 
calculate  pottucial  energy  curves  for  the  electronic  states  of  NF  has 
not  been  successful.  The  MNDO  program  produces  electronic  states  which 
cannot  be  identified  with  the  proper  diatomic  syrrmetry  for  certain  values 
of  the  intemuclear  coordinate.  Also,  one  would  like  to  include  many 
more  configurations  in  the  configuration  interaction  than  the  program 
allows.  Our  recommendation  is  that  potential  surfaces  for  small  molecules 
be  left  to  the  more  sophisticated  ab  initio  schemes. 

Heats  of  formation  and  geometry  of  boron-containing  species  of 
rocket  plume  interest  with  MINDO/3  are  remarkably  good  in  sane  cases  and 
rather  bad  in  others.  The  MNDO  program  results  are  needed  in  order  to 
make  an  assessment  of  the  usefulness  of  this  technique  for  these  appli- 
cations. 

Calculations  of  kinetic  parameters  for  TOT  decomposition  have  really 
only  begun.  The  initial  attempts  have  already  produced  an  estimate  for 
a first-step  activation  energy  value  for  the  rate-determining  step. 

These  calculations  should  prove  extremely  useful  in  deciding  among 
potential  mechanisms  for  TOT  decomposition. 

Initial  attempts  to  correlate  oxidation  potentials  for  aromatic 
molecules  with  the  calculated  molecular  orbital  energies  have  been  very 
successful.  These  correlations  can  hopefully  be  used  to  identify 
anomalous  systems  in  which  appreciable  solute-solvent  interactions  are 
occurring. 
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